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Классификация радиологических 
изображений в пространстве 
признаков ядер Вольтерра 


Рассмотрены информационные возможности нового метода классификации радиологических изображений. 
Суть метода основана на разделении анализируемых компонент на две составляющие, при этом 
высокочастотная компонента рассматривается в качестве входной, а низкочастотная -— в качестве 
выходной составляющей нелинейного формирующего фильтра, описываемого ядрами Вольтерра 
второго порядка. Представлены результаты реальной проверки работоспособности метода. 


Введение 


Одной из важных задач медицины является ранняя диагностика и профилактика 
остеопороза, приводящего к увеличению пористости костной структуры из-за «вымыва- 
ния» кальция вследствие нарушения поцессов обмена веществ в организме. В настоящий 
момент диагностика возможна только на основе анализа радиологических изображе- 
ний (рентгенограммы, компьютерная томография). Однако информационные возможно- 
сти визуального анализа радиологических изображений серьезно затруднены влиянием 
трех основных факторов, потенциально приводящих к многочисленным ошибкам диаг- 
ностики: 

1) недостаточной чувствительностью рентгеновских лучей к небольшим локаль- 
ным вариациям изменений плотности костных структур; 

2) недостаточной чувствительностью человеческого глаза, воспринимающего, в 
соответствии с законом Вебера [1], лишь 2% перепады яркости, причем это значение 
существенно зависит от окружающего фона; 

3) недостаточной пространственной разрешающей способностью радиологических 
изображений (в рентгеновской томографии решение этой проблемы связано с необхо- 
димостью увеличения числа проекций, однозначно приводящих к увеличению степени 
радиологической нагрузки на организм, а в простых рентгеновских аппаратах это свя- 
зано с проблемой фокусировки рентгеновских лучей [2]). 

Все вместе взятое это означает, что снятые, например, с промежутком в один год 
радиологические изображения могут фактически не различаться с точки зрения ви- 
зуального анализа, тогда как на самом деле костная структура может претерпевать су- 
щественные физико-биологические изменения. 

Целью данной работы является демонстрация информационных возможностей 
принципиально нового подхода к задаче классификации радиологических изображений 
костных структур, основанного на использовании признаков ядер Вольтерра, применяю- 
щихся для анализа нелинейных систем. 
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Математическая модель костной структуры 


Учитывая отмеченные проблемы прямого визуального анализа и классификации 
радиологических изображений костных структур, пораженных остеопорозом, «тонкую 
диагностику», по нашему мнению, целесообразно проводить на основании анализа осо- 
бенностей яркостно-пространственных срезов изображений костных структур. 

Базируясь на общих физических принципах, томографическое изображение кост- 
ной структуры /(х, у) можно представить в виде 


Их, у) = 5(%, у) + 2(х, у) +п(х, У), (0 
где 5(х, у) представляет собой низкочастотную нестационарную компоненту, харак- 
теризующую изменение общего фона изображения костной структуры; 2(х, у) — вы- 
сокочастотная стационарная составляющая, характеризующая динамику изменения 
пористости структуры, которую, собственно говоря, и можно рассматривать в виде 
аналога изображения текстуры; и(х, у) — измерительный шум. 

Нейтрализация влияния измерительного шума базируется на использовании стан- 
дартных методов фильтрации, поэтому основной является задача выделения компо- 
нент 5(х,у) и 2(х,у) из анализируемого изображения /(х, у). С этой целью нами 
был использован селективный метод автоморфного отображения [3], причем в качестве 
5(х, у) использовалось первое, а в качестве 2(х, у) — второе собственные изображения 
автоморфного отображения (рис. 1). 


Рисунок 1 — Томографическое изображение тазобедренного сустава (верхний ряд — 
здоровый сустав, нижний ряд — больной): а) — оригинал; 6) — 5(х, уж); в) -— 2(Х, Ур) 


На рис. 1 показаны томографические изображения здорового и больного тазобед- 
ренного суставов и яркостно-пространственные срезы 120-х строк, составляющих 5(х, у) 
и 2(х, у) соответственно. 

В соответствии с рассмотрением костной структуры в виде аналога пористой 
текстуры компоненту 2(х,у„., ) можно выразить в виде 


20) =2(х ,у.„)= У, бб -ху)+а(9, ©) 
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где Б(х) рассматривается в виде функциональной зависимости, характеризующей оги- 
бающую пространственного изменения яркости пористой структуры; х, — расстояние 


между порами; а(х) — шумовая непредсказуемая компонента. Поскольку в больной 
кости размер пор и их взаимное расстояние должны быть больше, чем для здоровой 
кости, то может показаться, что решение задачи диагностики возможно на основе 
использования методов спектрального анализа, поскольку для больной кости должна 
доминировать низкочастотная компонента, что, однако, не соответствует действитель- 
ности, как это следует из рассмотрения нормированных энергетических Фурье-спектров 
компонент 2(х) (рис. 1в), представленных на рис. 2, в предположении, что х, -х,/=1. 


1 


1 


0.9 0.9 
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0.6 |---- 0.6 --- 
0.5 0,5 


0.4 ------- 0.4 
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Рисунок 2 — Нормированные энергетические Фурье-спектры высокочастотных 
компонент, представленных на рис. 1в): а) — здоровый сустав; 6) — больной 


Из анализа рис. 2 следует, что интенсивность спектра |2 (Г ) на частотах / > 0,1 


для больной кости действительно в два раза больше, чем для здоровой, однако это 
отличие не носит качественного характера, поскольку максимумы обоих характерис- 
тик |2 (Г ) расположены на одинаковой частоте / = 0,07 ‚, что не позволяет класси- 
фицировать структуры по виду их спектральной характеристики. 

Исходя из вида спектральных характеристик |2 (Г) 


‚ модель костной структуры (2) 


можно представить в виде 


2(х)=2 + У, соз(2л/,х+ф,) + а(х), (3) 


К=1 
где 2, характеризует среднее значение яркостного фона, Б,, }, и ф, — амплитуды, частоты 
и фазы скрытых периодичностей (/ = 0,07 — доминирующая гармоника, отражающая 


среднюю «плотность» заполнения порами костной структуры). 


В свою очередь, модель (3) может быть представлена в виде модели авторегрес- 
сии [4] 
Р=2К 


2(х) = 2% + У`4(02(х —[)+а(>х), (4) 


где 4(Г) - вектор коэффициентов линейного предсказания, а Р’ -— порядок модели ли- 
нейного предсказания. 
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В рамках модели (4) предполагается, что а(х) является входной последователь- 
ностью линейного фильтра, формирующего выходную последовательность 2(х), при- 
чем под а(х) понимается ошибка прогноза, т.е. 


а(х) = 2() 209, (5) 
которая должна удовлетворять условиям 
Е{а(х)}=0; Е‘а(х,)а(х,)}= 026-10, (6) 


где ЕЁ — оператор математического ожидания. Такой подход к модели пористой струк- 

туры открывает возможность вычисления параметрической спектральной характерис- 
^ 

Е(Л) 


ТИКИ на основе выражения [4] 


|2 (Г) = ИНН ‚ (7) 
П-У'а@Фехр(- 2х) 


обеспечивающего более высокую разрешающую способность спектрального анализа 
по сравнению с классическими методами Фурье-спектроскопии. 

На рис. 3 представлены параметрические спектры для данных на рис. 1в) и по- 
рядке модели линейного предсказания Р =14. 
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высокочастотных 


2() 


компонент (рис. 1в): а) — здоровый сустав; 6) — больной 


Рисунок 3 — Нормированные параметрические спектры 


2() 


ляется доминирующая гармоника ы ‚ а сами спектры не имеют каких-либо качест- 


Из рассмотрения рис. 3 следует, что на зависимости ‚ как и на рис. 2, выде- 


венных отличий. 
Из общего анализа представленных результатов можно сделать следующие выводы. 
1. Среднее количество пор как для больного, так и для здорового участков оста- 
ется примерно одинаковым, что не позволяет использовать стандартные методы спек- 
трального анализа для диагностики потенциального развития остеопороза на ранней 
стадии заболевания. 
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2. Если исходить из того, что развитие остеопороза приводит к увеличению раз- 
меров пор, без изменения их среднего количества в пределах каждого выборочного 
яркостно-пространственного среза, то задачей диагностики является выявление про- 
странственного взаимодействия пор, что диктует необходимость перехода от линей- 
ных моделей типа (2) и (3) к классу нелинейных моделей. 


Классификация как задача нелинейной фильтрации 


Вольтерра показал [5], что для стационарной нелинейной системы с конечной 
памятью соотношение между входом и(Г) и выходом у(Г) можно записать как 


У(В = [1 (т)и(Е-т)ат + [ [№ (т. )и@-т )и@-т.)ат, ат, + 


2—8 


тт) (8) 


где й (т); И, (т.т); №, (т,т,,т.) — ядра Вольтерра, являющиеся симметричными 


функциями, причем член первого порядка представляет собой интеграл свертки ли- 
нейной системы. Функциональный смысл ядер Вольтерра заключается в проведении 
аналогии между ними и взаимодействием нескольких импульсов. Поскольку ядро пер- 
вого порядка рассматривается как модель, с помощью которой учитывается влияние 
прошлых значений сигнала на реакцию в настоящий момент, то ядра Вольтерра выс- 
ших порядков можно представить как модель, описывающую зависимость значения 
выходного сигнала в настоящий момент от взаимодействия прошлых значений вход- 
ного сигнала. 

Алгоритмы вычисления ядер Вольтерра подробно описаны в [6]. Однако при- 
менительно к рассматриваемой в работе задаче классификации костных структур, 
применение модели (8) наталкивается на одну трудность — что рассматривать в ка- 
честве входного сигнала и(Г) =и(х)? В авторегрессионной модели (4) сигналом на 
входе формирующего фильтра рассматривается ошибка прогноза а(х), однако для мо- 
дели (8) такой подход неприемлем. Поэтому решение задачи базируется на использо- 
вании двух гипотез: 1) высокочастотная компонента 2(х, у) рассматривается в виде 
инновационной составляющей на входе формирующего фильтра для низкочастотной 
компоненты 5(х, у); 2) формирующий фильтр для компоненты 5(х,у) является нели- 
нейным и описывается системой Вольтерра второго порядка. 

Поскольку процедура «прямого» спектрального анализа зависимостей 2(х)не 
приводит к каким-либо качественно различающимся результатам, то основная идея 
предлагаемого подхода заключается в анализе именно ядер Вольтерра второго 
порядка й(х,,х,). 

На рис. 4 представлены результаты, полученные для данных, показанных на рис. 1. 

Если рассматривать изображения ядер Вольтерра второго порядка в виде 
изображений «виртуальных аналогов» текстур, то из анализа рис. 4а) следует, что 
здоровая кость является мелкопористой, а больная — крупнопористой, что, собственно 
говоря, и соответствует модели развития остеопороза. С физической точки зрения это 
означает, что здоровая кость имеет меньшую «пространственную память» Дх = х, —х,, 
по сравнению с больной, что следует из анализа диагональных срезов ядер Вольтерра 
на рис. 4 6) и проявляющихся в наличии интенсивной высокочастотной компоненты у 
здоровой кости (рис. 4 в), тогда как у больной кости доминируют низкочастотные состав- 
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ляющие. Таким образом, если прямой спектральный анализ яркостно-пространствен- 
ных срезов не дает положительных результатов, то Фурье-анализ диагональных срезов 
ядер второго порядка выявляет «тонкое» различие в структурах больной и здоровой кости 
и может служить количественным критерием этого отличия (у здоровой кости частота вто- 
рой доминирующей гармоники /, = 0,17 ‚ ау больной — Г, = 0,06). 


09-1 -- 
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Я И Я И | И | | | 0 + - | 
002 4 0 Пл 98 10 120 140 160 180 20 0 0.05 04 0.45 02 0.25 03 035 0.4 0.45 0.5 


отсчеты частота 


|] О А 
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х, отсчеты частота 
а) 6) в) 
Рисунок 4 — Результаты классификации радиологических изображений, 
представленных на рис. 1 (верхний ряд — здоровый сустав; нижний — больной): 
а) й(х,,х,); 6) диагональный яркостно-пространственный срез «а»; 
в) энергетический Фурье-спектр диагональных срезов 


Заключение 


1. Метод анализа радиологических изображений в пространстве признаков ядер 
Вольтерра позволяет проводить качественную и количественную классификацию кост- 
ных структур на предмет выявления заболевания остеопорозом. 

2. Теоретическая особенность метода заключается в рассмотрении анализируемого 
изображения в виде суммы двух компонент: низкочастотной — характеризующей из- 
менение общего фона костной структуры, и высокочастотной — характеризующей 
динамику изменения пористости структуры, причем эта компонента рассматривается 
в качестве входа нелинейного формирующего фильтра, описываемого системой Воль- 
терра второго порядка, для низкочастотной фоновой составляющей. 
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3. Качественная диагностика возможна на основе анализа изображений ядер вто- 
рого порядка при рассмотрении их в качестве виртуального аналога изображений 
текстур, а количественная — на основе спектрального Фурье-анализа их диагональ- 
ных срезов. 

4. Проведенные исследования показывают, что предложенный подход имеет зна- 
чительный потенциал дальнейшего развития применительно к задачам классифика- 
ции, например, на основе дополнительного применения методов биспектрального 
анализа. 


Литература 


1. Ргай \\.К. Плеиа! Ппасе Ргосеззте: РК$ ше, Тьиа Еф@оп / Ргай \У.К. . — М.У : оба УШеу & 
501$, 2001. — 783 р. 

2. Физика визуализации изображений в медицине / под ред. С. Уэбба. - М. : Мир, 1991. -— Т. 1. - 
406 с. 

3. Ахметшина Л.Г. Адаптивная фильтрация шумов в сигналах и изображениях: метод селективного 
сингулярного разложения автоморфного отображения / Л.Г. Ахметшина // Искусственный интеллект. — 
2005. — №3. -— С. 328-335. 

4. Марпл С.Л. Цифровой спектральный анализ и его приложения / Марпл С.Л. - М. : Мир, 1990. — 
584 с. 

5. Виззсапс /.7. Апа]уз1з о попПпеаг зу%ет ул шире триб / 7.7. Виззеапе, Г.. Ебштап, 1.М/. Отават // 
Ргос. ТЕЕЕ. - 1974. — Уо/. 62. -Р. 1088-1119. 

6. Мармарелис П. Анализ физиологических систем: метод белого шума / П. Мармарелис, В. Мар- 
марелис . — М. : Мир, 1981. —480 с. 


О.М. Ахметшин, К.О. Ахметшин 

Класифжкащя радюлойчних зображень у простор! ознак ядер Вольтерра 

Розглянуто 1иформащйн! можливост! нового методу класифжацй радологичних зображень. Суть 
методу полягае у розподл! аналзованих компонент на дв! складов1, при цьому високочастотна 
компонента розглядаеться як вх1дна, а низькочастотна — як вихдна компонента нелшййного 
формуючого фальтра Вольтерра другого порядку. Представлен! результати реально! перевфки 
працездатност! методу. 
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Са$$ШсаНоп о! Ва4101001са! Ппасе$ ш а 5расе оЁ УоНегга Кегпе! Ееафиге$ 

шЮгпайоп роззоШНез оРа пе\у тео4 с1аз$1сайоп оРтаФоос1са| ппазэез аге сопз14егед. Тве тат 14еа 
оЁ Фе те@фод 1$ Базе оп зерагайоп оЁ апа]уте итазе оп мо сотропет: №1е6-Недаепсу сотропеп! 15 
три! апа 1о\/-Яедиепсу опе 1$ ори попПпеаг Юглте УоЦегга НШег оЁ зесоп4 ог4ег. ТБе гезиЙ$ оР геа] 
4е5 Ня оРфе тефо4 роз Швез аге ргезеще4. 


Статья поступила в редакцию 02.07.2009. 
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